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THE MOMENT OF INSTABILITY FOR 
INTERNAL SOLITARY WAVES 

ANDREAS KLAIBER 


Abstract. In this note, we define a moment of instability m(c) for 
internal solitary waves in continuously stratified fluids, which seems not 
to have been done before. To underline the suitability of the proposed 
m(c), we identify the relation m"(c) = 0 as a formal Fredholm condition, 
and we show that m"(c) displays a definite sign for small-amplitude 
waves. 


1. Introduction 

Internal solitary waves (ISWs) are ecologically important since they are 
involved in mixing mechanisms and energy transport in lakes and oceans 
In this context, a widely used mathematical model consists of 
the 2D Euler equations for incompressible, inviscid fluids with non-constant 
density. This model is given by the equations 

(1.1a) p t = -up x -vp y , 

(1.1b) ut = — uu x — vu y — —, 

D 

(l.lc) v t = -uv x - vv y - -J - g, 

complemented by the incompressibility constraint 
(l.ld) 0 = u x + v y , 

the boundary conditions 

(l.le) v(t,x, 0) = 0 and v(t, x, 1) = 0, 

and the far-held conditions 

(I lf) {p,u,v,p)(t,±oo,y) = (p(y),0,0,p(y)), 0 < y < 1. 

In (11.111 . density p, velocity (u,v), and pressure p are functions of time t. 
horizontal position x € R and vertical position y £ [0,1], and the constant g 
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denotes acceleration due to gravity. The far-held (p(y),0,0,p(y)), itself an 
x- and f-independent solution of (ll.lall - dl.lcl) with 

p : [0,1] -A- (0, oo) differentiable with p'{y) <0, 0 < y < 1, 

ry 

and p(y) = ~g p(rj ) d?y, 

Jo 

is called the quiescent state. Travelling wave solutions 

(p, u, v,p)(t , x, y) = (p, u, v,p)(x — ct, y), with some c > 0, 

of (11.11) are called internal solitary waves (ISWs) of speed c; we refer to 
[smaii] for mathematical results on their existence. 

In order to study the stability of ISWs, one could start from the lineariza¬ 
tion of (11.11) about a given ISW, as done by the author in [TUI HI] to find an 
Evans-function approach to stability. 

A different general approach to investigate the stability of solitary waves 
is based on the moment of instability, see e. g. [6], and references therein. 
Here, we want to establish the moment-of-instability (MOI) route to stability 
of ISWs. 


2. Definition of a moment of instability m(c) for ISWs 

According to [3], the Euler equations CEO for stratified fluids possess 
a Hamiltonian formulation. In terms of the density p, the vorticity-like 
quantity a, and the associated streamfunction ip, defined as the solution of 


cr = -V ■ (pVip), with ip | y=0i i = 0, 

this can be formulated as 

V 


( 2 . 1 ) 


a 


dt\'\ = J(p, a) (H - cl) (p, a) 


in the co-moving frame t, x = x — ct, y with writing x instead of x, where 
the Hamiltonian Ji — cl is composed of the energy functional 

(2.2) H{p,a)= [ [ \p\S7ip\ 2 + gy(p-p)dydx 


and the momentum functional 
(2.3) X{p,a) = 


ya d y dx, 


and J = J(p,a) denotes the (state-dependent!) skew-symmetric operator 
(2.4) J(p,c j) = 


0 Pxdy T Pydx 

Px&y T Py(Jx T xIJy 4“ CTyd^ 

The Hamiltonian formulation (12.11) . however, does not directly yield a varia¬ 
tional principle due to the non-invertibility of J. Concretely, as a stationary 
solution of (12.11) an ISW ( p c , a c ) satisfies 


( 2 . 5 ) 


0 = J(p c ,a c ) (H-cX) ( p c ,a c ) 






THE MOI FOR ISWS 


3 


but, as a little calculation reveals (see, e. g., m p- 35 ]), 

( 2 . 6 ) (n - cl) (p c , a c ) = (v y - |i v ^ c|2 ) + 0 , 

i. e., (p c , a c ) is not a critical point of H — cl\ 

This issue can be overcome by modifying H — cl without spoiling the 
Hamiltonian structure. In fact, taking the quantities 

(2.7a) AH(p, a) := — f [ g\[ p~ l (g) dg\a dy dx, 

JW JO ^ J p(y) J 

(2.7b) Al(p,a):=- [ [ p~ 1 (p)adydx , 

Jr Jo 

it is easily verified that 

(2.8) J(p, a) (A H - cAZ)' (,p , a) = 0 and (U - cl)' (p c , a c ) = 0 
with % :=% + AH and 1 = 1+ A1. Therefore, replacing H — cl with 


U-cl = 


ff 


-pIVY^ + W {y - p 1 (g)} dgdydx 


' p{y) 


— c 


{ y-p 1 (p)}°dydx, 


results in a modified Hamiltonian formulation such that ISWs are, indeed, 
critical points of the Hamiltonian. This was already noticed by pH I] but, 
as far as the author is aware, has not been used in connection with the 
stability of ISWs. For background material on so-called Casimir functionals, 
for which AH — cAl is an example, their systematic derivation and their 
use in hydrodynamic contexts, see |T| and references therein. 

Now, we are in a position to define the moment of instability for ISWs in 
the usual way: 


(2.9) m(c) := (H — cl) ( p c ,cr c ). 

Since (H — cl) (p c , er c ) = 0 by construction, we immediately have the usual 
relation 

m"(c) = ^ (H - cl) (p c , u c ) = -^X(p c , u c ). 

In the rest of the paper, we study this m(c). In Sec. 2 we show that 
m"(c) < 0 for ISWs of sufficiently small amplitude. In Sec. 3 we show in a 
quite general situation, which covers ours, that the condition m"(c) = 0 can 
be read as a formal Fredholm condition. 

Acknowledgment. I thank Heinrich Freistiihler for drawing my attention to 
the fact that the meaning of m"(c) = 0 can be understood without reference 
to a possibly existing Evans function, and for stimulating discussions on this 
topic. 
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3. Proving m"(c) < 0 for small ISWs 
For small wave^l, we have mm 
c = cq + £ 2 , 

^ c (x, y) = £ 2 A{£x)yo{y) + 0(e 4 ), 

P c {x,y) = p{y) - — e 2 A(ex)p'(y)p 0 (y) + 0(e 4 ), 

co 

where 

A"(X) = --A(X ) - -A(X) 2 and (p(y)^o(y))' = ^p\y)p 0 {y). 

S S Cg 

With these expressions at hand, it is straightforward to evaluate m"(c). 


l{p c 1 a c ) = - 
c 

1 

c 


p c |V#f dydx 


1 


p(y) -c A(ex)p (y)ipo(y) + 0(£ ) 

co 


x ((£ 3 A'(£x)ip 0 (y)) 2 + ( £ 2 A(£x)ip' 0 (y)) 2 + 0(e 5 )) dydx 

— [ [ p(y)M£x) 2 <Po(y) 2 dydx + 0(£ 5 ) 

Co JR Jo 

[ A(£xfdx [ p(y)<p' 0 {y) 2 dy + 0(e 5 ) 

Jr Jo 

— [ A(X) 2 dX f p(y)p' 0 {y) 2 dy + 0(£ 5 ) 

c 0 JR JO 

= A" (c — c 0 )^ + O ^(c - c 0 )^ 


with the finite, positive constant 
K := - 

Hence, we derive that 
d 


co Jr 


A{X) 2 dA [ p(y)p' 0 (y) 2 dy > 0. 
Jo 


m"{c) = --^l[(p c ,a c )\ = --K (c- c 0 ) 2 + O ((c- c 0 )2^ < 0 
holds for 0 < c — co 1, i. e., for sufficiently small waves. 


4. Characterizing m"(c) = 0 as a Fredholm condition 

To simplify the notation, we write 4> = (p C i (jC ) f° r the ISW in the follow¬ 
ing. In the situation above, differentiating the profile equation 

(4.1) (A - cX)'(^) = 0 
with respect to the position yields 

(4.2) ( H - cZ )" ( 4 )^ = 0 , 

* We assume here the genericity condition J 1 p(y)ifio(y) dy 7^ 0 which is necessary for 
the validity of the approximate expressions; cf. [5]. 
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while differentiating it with respect to the speed results in 

(4.3) 

Eqs. (14.21) . (14.31) give 

(4.4) JW~cZ)"(<j>)^ = 0 
and 

(4.5) = 

As 

(4.6) Jl'W = 

eqs. (14.4p and (14.51) state that 0 is an at least double eigenvalue for 

u = J{U - cl)'(u), 

with ^ as an eigenfunction and ^ as a first-order generalized eigenfunction. 
Now, a second-order generalized eigenfunction ■*/> would solve 

(4.7) = 

According to the Fredholm alternative, eq. flED has a non-trivial solution 
if and only if its right hand side ^ is orthogonal to the solution % °f the 
adjoint homogeneous equation 

(4.8) 0 = (J(H - cl)"m*x = -m - cl)"^)J) X . 

As (14.21) and (14.61) imply 

(4.9) 0 = —((W. — cl)"((j>)J)l'((f>) and thus x = I / (</>), 
the existence of ijj ^ 0 consequently is equivalent to 

(4 - 10) 0 = l X ^) = ( Xm f), 

i. e., vanishing of the moment of instability. 

Remarks, (i) The above argument slightly varies the one given by Zum- 
brun in [16] (Sec. 1, between the statements of Corollary 1.4 and Remark 

1-5). 

(ii) This argument literally applies to the situation of Grillakis et. al. P3 b y 
changing to their notation 

E = 'H,Q=l,J = J,(j)= ( p c , a c ),T'{ 0) = d x . 


Hence, it can be applied to various contexts that fall into this class. 
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